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ABSTRACT 

The upcoming generation of cosmic microwave background (CMB) experiments face a 
major challenge in detecting the weak cosmic B-vnode signature predicted as a product 
of primordial gravitational waves. To achieve the required sensitivity these experiments 
must have impressive control of systematic effects and detailed understanding of the 
foreground emission that will influence the signal. In this paper, we present templates 
of the intensity and polarisation of emission from one of the main Galactic foregrounds, 
interstellar dust. These are produced using a model which includes a 3D description of 
the Galactic magnetic field, examining both large and small scales. We also include in 
the model the details of the dust density, grain alignment and the intrinsic polarisation 
of the emission from an individual grain. We present here Stokes parameter template 
maps at 150 GHz and provide an on-line repositor3|t]for these and additional maps at 
frequencies that will be targeted by upcoming experiments such as EBEX, Spider 
and SPTpol. 

Key words: cosmic microwave background, polarisation experiments, foregrounds, 
5-modes, gravity waves 



1 INTRODUCTION 

The next round of cosmic microwave background (CMB) 
experiments are all targeting measurements of CMB polari- 
sation. The CMB polarisation field can be decomposed into 
curl-free ^-modes and curl-like B-modes (see for example 
[Kamionkowski, K osowsky fc Stebbins' (1997)), however to 
date only ^-modes have been observed ([Kovac et al. 2002 



Readhead et al.|2004|[Sievers et al.|2007|[Barkats et al 



Bischof f et al.|20 08"Leitc h et al.|2005||Montroy et al. 



[Page et al.| 



2007 



Piacentini et al., 2006; Wu et al.||2007 
Brown et al.||2009| [Chiang et al.||20Tor 
tion||2010). Current work centres on the search for the tiny 



2005 



2006 



frequencies being targeted by EBEX, Spider and SPTpol. 
Galactic foreground emission is expected to contribute sig- 
nificantly to the polarised microwave emission across the sky, 
and may well dominate the CMB gravitational-wave signal 
at all frequencies. In order to achieve their scientific goals, 
forthcoming CMB polarisation experiments require in-depth 
knowledge of this polarised Galactic foreground emission. 
The Planck mission ("The Planck Collaboration||2006| will 



QUIET Collabora-I 



amplitude S-modes since a key prediction of inflation is the 
generation of a stochastic background of gravitational waves, 
which on large angular scales are the only contribution to 
the CMB B-uiode component. Experiments involved in this 



provide maps of the polarisation of interstellar dust, allow- 
ing tests of the structure of the Galactic magnetic field. It 
should also provide an insight into grain alignment mecha- 
nisms. 

Of utmost importance will be the accurate separation of 
foreground emission from the CMB signal. Component sep- 
aration has been considered in, for example, Brandt et al 



Spider (Filippini et al. 2010), SPTpol (McMahon et al. 



et al.|2010|), AC Tpol ( [Niemack et al.|2010|, POLAR (Keat- 
ing et al.||19'98t, PQLARBEAR jThe Polarbear Collabora- 



tion et al. 2010) and BRAIN (Charlassier & the BRAIN 



search include EBEX (|Reichborn-Kjennerud et al. 2010), et al.| (2009). In order to test and assess methods for sep- 



20091, PIP ER (Laze ar et al.|2011t, ABS (lEs singer-Hileman 



(1994}; Eriksen et al. (2006); Kogut et al. (2007); Stompor 



Collaboration 2008). We focus on providing templates at 
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arating out the contribution realistic foreground templates 
will be required. Since however, few data exist at this time 
at the observing frequencies in which the upcoming exper- 
iments are operating, one must resort to modeling of fore- 
ground emission by extrapolating the information from ex- 
isting data. Furthermore, experiments which observe por- 
tions of the sky close to the Galactic plane will find the fore- 
ground emission is very bright in comparison to the CMB. 
The presence of this bright emission in the data may affect 
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the performance of the observation and the analysis strategy 
an experiment uses. Another important role of foreground 
modeling, therefore, is informing the planning and proposal 
stage of any experiment. 

Unfortunately, these foregrounds are poorly constrained 
by current data and poorly understood, particularly above 
around 90 GHz, where the CMB emission is strongest and 
where many new CMB experiments will operate. At these 
frequencies, the foreground emission is expected to be dom- 
inated by thermal emission from interstellar dust. A review 
of the basic physical processes whereby aligned dust grains 



generate polarisation is given in Lazarian & Cho (2003) 
Such emission is known to be polarised both through di- 
rect measurements (Ponthieu et al.|2005 Kogut et al.|2007 



Benoit et al. 2004 Bierman et al.|2011) and through obser- 



vations of the polarisation of starlight ( Heiles|1996 Fosalba 
eraL][2002t . 

This polarisation arises due to the presence of a mag- 
netic field in the Galaxy. The dust grains are generally non- 
spherical, and preferentially emit radiation polarised along 
their longest axis. Mechanisms exist which align these grains 
with this axis perpendicular to the Galactic magnetic field, 
leading to net linear polarisation. 

Polarised foregrounds also include polarised emission 
from synchrotron. Synchrotron emission, generated by the 
gyration of cosmic ray electrons in the Galactic magnetic 
field, is intrinsically polarised and constitutes the main po- 



larised foreground at lower frequencies (Page et al. 2007). 



However, emission from thermal dust dominates synchrotron 
emission at the frequencies considered in this paper and 
therefore we concentrate on modeling emission from only 
thermal dust. Foreground radiation also includes free-free 
and spinning dust emission, however we assume both of 
these signals to be unpolarised and so do not consider them 
further here. Evidence for this has emerged recently, with 



Macellari et al. (2011 ) showing that free- free emission is un- 
polarised, setting an upper limit on the free-free polarisation 
fraction of 3.4% at the 2cr level. They also show that spin- 
ning dust emission has a low polarisation fraction. Upper 
limits on the polarisation fraction of spinning dust emission 
in molecular clouds have been obtained by Dickinson, Peel| 
k Vidal| ( [20T r) and Lopez- Carabal lo et~al] ( 2011| ), who find 
low levels of polarisation. If there is a similar level at higher 
Galactic latitude, then this foreground is unimportant in 
terms of component separation. 

The model in this paper was introduced by |0'Dea| 



( 20091 and first applied in |0'Dea et al.| ( |2011| ) for the pur- 
pose of studying the impact of polarised foregrounds on Spi- 
der's ability to detect B-mode polarisation. Here we give 
a more detailed explanation of the model and present a 
number of full-sky template maps at various frequencies. 
Arch eops ( [Benoit et al.|2004| and BICEP ( [Bierman et al. 



2011) have made the highest signal-to-noise maps of the 



dust polarisation at low Galactic latitudes and have exam- 
ined the properties of the polarisation fraction and angle. 
These results cannot be relied upon to calibrate large scale 
models of the polarisation at higher latitudes since the po- 
larisation properties near the Galactic plane will depend on 
complex structure that is not included in models such as 
the one presented in this work. Here, both polarisation am- 
plitude and angle are modeled internally and our templates 
are scaled such that the polarisation fraction corresponds 



to a nominal value when averaged over the maps with the 
Galaxy masked out. The Archeops and BICEP maps have 
not been made public and a quantitative comparison of our 
templates with these observations is not possible. [Fauvet,| 
L. et al. ( 2011| have developed a similar model of both the 
polarisation of thermal dust and synchrotron radiation and 
compared with the WMAP K-band and Archeops 353GHz 
data, however they have not released any templates based 
on their model. 

This paper is organized as follows. We begin in Sec- 
tion |2]T] by describing the model of the dust total intensity. 
We focus on modeling the Galactic magnetic field (which is 
made up of a large-scale field from the Galactic disk and a 
small-scale field due to turbulence in the interstellar plasma) 
in Section [2^ We include in Section [23] a description of the 
large-scale dust density, as well as details of the treatment 
of dust grain alignment and the intrinsic polarisation of the 
emission from an individual grain. In Section [2.4[ we show 
how these fields are modeled in three-dimensions and are 
used to perform a line-of-sight integral to the centre of each 
pixel to form the final Stokes parameter maps. In Section |3] 
we discuss the effective resolution of the maps produced us- 
ing our model and the relevant physical scales introduced by 
various effects. In Section [4] we describe the final template 
maps and summarise our conclusions in Section |5] 



2 FOREGROUND DUST MODEL 

2.1 Dust Total Intensity 

Although few data are available regarding the polarised 
emission from dust, the same is not true of its total inten- 
sity. In particular, the IRAS satellite observed this emission 
across the sky at lOO/xm and 240/im, close to the peak in the 
dust emission. By constraining physically-motivated extrap- 
olations of these observations using further data, [Finkbeiner'^ 
Davis & Schlegel (1999[ hereafter FDS) provided models of 
the emission at microwave wavelengths. At 94 GHz, these 
models have been shown to agree well in terms of morphol- 
ogy with the WMAP observations with some minor struc- 



tural differences on the Galactic Plane (Gold et al. 2009) 



However, there are indications that in terms of amplitude 
the WMAP dust template fit coefficients differ by about 



30%. Bennett et al._ (2003 ) suggest that this is possibly due 
to the degeneracy that exists between the strongly corre- 
lated dust and synchrotron emission components in the si- 
multaneous fit of their externally derived template maps to 
WMAP data. 

In the higher frequency bands relevant to experiments 
observing above ^ 90 GHz, data are more limited but agree 



well with the FDS predictions (Culverhouse et al. 2010 



Veneziani et al. 2010). We will use this model (to be pre- 



cise model number eight in FDS) to trace the total intensity 
of the dust emission. We exploit the full, 6.1 arcminutes, 
resolution of the IRAS data by pixelising the dust intensity 
on HEALPIX (Gorski et al.[[2005) maps of A^side=1024. 



2.2 Galactic Magnetic Field 

The degree and direction of polarisation of the dust emis- 
sion are highly dependent on the Galactic magnetic field. As 
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Figure 1. Cartesian projection of the BSS (left panel) and LSA (right panel) large-scale magnetic field models, showing a slice through 
the Galactic plane observed along the positive z-axis. The filled circle represents the position of the Sun in relation to the Galactic 
centre. The alignment and magnitude of the magnetic field are shown as headed ticks with the Bq = 2 fiG scale represented at the top 
of each panel. The red solid lines are the density contours, in steps of 0.1, of the dust density model nd(r,z) (see Section [z!3| used in 
the line-of-sight integration. The dust density is normalised to 1 at the Galactic centre. 



the observed polarisation is the sum of many independent 
regions along the line-of-sight, it is sensitive to the three- 
dimensional structure of this magnetic field. Therefore, to 
proceed we first consider a three-dimensional model of the 
Galactic magnetic field (and the other necessary Galactic 
constituents and physics) and then set the polarisation de- 
gree and direction through the appropriate line-of-sight in- 
tegrals. 

Away from the Galactic center, the Galactic magnetic 
field is usually considered to have two near-independent 
components: a large-scale coherent field associated with the 
Galactic disk, and a small-scale field arising from turbulence 
in the interstellar plasma sourced by astrophysical events 
such as supernovae and stellar winds. The most informa- 
tive probes of these fields are Faraday rotation measures of 
pulsars and extra- Galactic radio sources f Haver korn et al. 
|2006[ I Han et al. 2006). Whilst there is general agreement 
that the large-scale field follows a spiral pattern, its detailed 
structure is still uncertain. 

When considering areas of sky at high Galactic lati- 
tudes, this uncertainty is unimportant as the dust is con- 
centrated in a thin disk about the Galactic plane, and so we 
only see emission within around 1 kpc or so of the Sun, a 
region in which the large-scale field is reasonably well char- 
acterised. However, experiments which will target a large 
fraction of the sky, possibly including part of the Galactic 
plane, will require a model of the large-scale field structure 
in the plane. 

Attempts have been made to constrain the proper- 
ties of the magnetic field using CMB polarisation measure- 
ments. jJafe^^T^T (2010) use an Monte Carlo Markov Chain 
(MCMC) approach to test components of a 2D Galactic 
field model using rotation measures and WMAP data in the 
plane of the Galaxy. Jansson et al. (2009) use rotation mea- 
sures and WMAP 5 year data to fit for parameters in com- 
mon 3D models for the Galactic magnetic field. We choose 



two of the most popular forms for the Galactic magnetic 
field and provide templates for both of these models. 



2.2.1 Large-Scale Magnetic Field 

One popular candidate is the Bi-Symmetric Spiral (BSS) 
(Han fc Qiao|1994) Sun et al.|2008) which can be written as 



Bp = —Bo cos \ ^ -\- thln^- I smpcos y , 



B^ = —Bo cos $ + V^ln — cos p cos X , 
Bz = Bo sin X ■ 



(1) 



Here, yo, $ and z are Galacto-centric cylindrical co-ordinates 
with $, the cylindrical longitude, measured from the di- 
rection of the Sun, p is the pitch angle of the field, ip — 
1/tanp, po defines the radial scale of the spiral, x — 
Xo tanh(2;/zo) parametrizes the amplitude of the z compo- 
nent and zo — 1 kpc. We use the parameters constrained 
in [Miville-Deschenes et aL (2008|: p — —8.5 degrees, po — 
11 kpc and xo = 8 degrees, with the field amplitude set to 
Bo — 3/iG, and take the distance between the Sun and the 
Galactic center to be 8 kpc. A diagram of the magnetic field 
orientation and magnitude in the BSS model is shown in the 
left panel of Figure [l] 

A number of other magnetic field models have been pro- 
posed in the literature. For comparison we also include the 
Logarithmic Spiral Arm (LSA) model introduced by Page] 
et al.l (|2007|) for use in cleaning of the WMAP data. The 
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model is defined as 



Bp = -Bq sin V^o + ^1 In 






cosx, 



B^ — —Bo cos ( '00 + '01 In J cos x 5 

B.^Bosinx, (2) 

with parameters obtained by fits to the WMAP K-band field 
directions; ipo = 27 degrees, ipi = 0.9 degrees, and x defined 
as in the BSS model but with xo = 25 degrees. The radial 
scale is also different in this model with pw — 8 kpc whereas 
the scale height is the same as above with zq — 1 kpc. There 
is no azimuthal dependence in this model. The right panel 
of Figure [l] shows a slice through the Galactic plane for the 
LSA model. 

Although both fields are unlikely to provide a full de- 



scription of our Galaxy (Men, Ferriere & Han 2008| Sun 



et al.|2008), they are sufficient for our current purpose as we 



do not require a precise template of the sky, only a reason- 
able approximation against which to test foreground separa- 
tion techniques and the performance of experiments in the 
presence of systematic effects. 

Both magnetic field models assume the field strength 
Bq is constant although there is weak evidence for some 
radial dependence (Han e t al.||2006 ). Any such dependence 
will not affect the polarisation model significantly and the 
overall radial dependence of the signal is determined by the 
exponential drop-off in the dust density which modulates the 
integrand along the line-of-sight. Field reversals may also 
be present in the spiral arms but, if sharp enough, will not 
contribute to the signal significantly. 



spectrum one finds that the power spectrum of a turbu- 
lent field is V{k) oc ^-(2+3iVd)/3 ^j^g^.^ ^^ jg ^j^^ number of 

spatial dimensions of the realisation. In this work, we gen- 
erate realisations of this power spectrum in order to model 
the three-dimensional magnetic field in real space by a Fast 
Fourier Transform (FFT). 

It is numerically intractable to generate a realisation of 
this turbulent field in three dimensions at sufficiently high 
resolution and to accommodate the entire sky, hence we re- 
sort to independent one-dimensional realisations along the 
line-of-sight to each pixel. This model ignores correlations 
across the sky, but properly incorporates the line-of-sight de- 
polarisation. We choose an injection scale of 100 pc, assume 
the dissipation scale is small and use the one-dimensional 
Kolmogorov energy spectral index of —5/3. 

For smaller patches of the sky, relevant for ground-based 
observations, a full, three dimensional realisation is feasible 
together with a higher angular resolution in the line-of-sight 
integrals. 



2.3 Dust Properties 

We model the large-scale spatial distribution of the dust 
density, n^, using a simplification of the model constrained 

in '. 



Drimmel fc Spergel (2001) 



rid = ^0 exp 



Pd 



sech' 



Zd 



(3) 



2.2.2 Small-Scale Galactic Magnetic Field 

The turbulent field is somewhat less well understood. When 
constraining the above large-scale field, M iville-Deschenesl 
|et al._ (2008 ) simultaneously fit a small-scale field with best- 
fit r.m.s amplitude Br.m.s. = 1.7/iG. Several different studies 
agree that the r.m.s. amplitude is similar to the amplitude of 
the large-scale field in the Solar vicinity ([Fosalba et al.|2()02 



Han et al. 2006), and so here we set B^. 



2/iG. 



Minter 



fc Spangler^ (^1996 ) examined the rotation measures of extra- 



Galactic sources across a small patch of sky and concluded 
that the data were consistent with Kolmogorov turbulence 
on scales smaller than 4 pc, assuming a statistically isotropic, 
homogeneous Gaussian field. On larger scales they found a 
somewhat flatter energy spectrum with an outer scale of up 
to 96 pc. Kolmogorov-type spectra up to kilo-parsec scales in 
the interstellar magnetic field and other interstellar plasma 



components have also been reported by other studies ( Arm- 



strong, Rickett fc Spangler|1995|[Lazarian fc Pogosyan 2000| 



Cho fc Lazarian||2008t ." 

Kolmogorov turbulence describes the energy distribu- 
tion among vortices of different size, with the amplitude of 
the turbulence related to the energy density at that posi- 
tion. Turbulent flow can be viewed as an energy cascade 
from larger to smaller eddies. At small enough length scales, 
known as the Kolmogorov length scale, energy is dissipated 
through viscous dissipation. A Kolmogorov spectrum is pro- 
portional to the rate of energy dissipation and the mag- 
nitude of the wavevector k. Using the Kolmogorov energy 



For consistency with the WMAP polarisation analysis 
( Page et al.|2007 ), we take the scale height Zd — 200 pc and 
the scale radius pd — "^ kpc. We do not attempt to model 
the small-scale variations in the dust density and tempera- 
ture here, which may also affect the polarisation degree and 
direction. Small-scale variations in the total intensity are 
included via the FDS model. 

The model also requires a description of the physics of 
grain alignment and of the intrinsic polarisation of the emis- 
sion from an individual grain. In general these are complex 
functions of the magnetic field and various properties of the 
grains. Recently, good progress has been made in describing 
the details of the alignment using the theory of radiative 
torques ( [Lazar ian & Hoang 2007| |Hoang fc Lazarian 2008 ) . 
However, it is still difficult to produce a well-constrained 



quantitative description to apply to our model ( [Lazarian fc 
Hoang|2009t . 



Instead, we describe the alignment in an integrated 
manner, without recourse to the details of a particular phys- 
ical mechanism. We assume that the polarisation direction is 
always perpendicular to the component of the magnetic field 
in the plane of the sky, and that the degree of polarisation 
depends quadratically on the magnetic field strength. This 



is similar to the behaviour assumed in Page et al. (2007) 



We follow this approach in providing our templates and do 
not attempt to account for any possible misalignment of the 
axis of orientation of the dust grains with the magnetic field 
lines, as is done in other work, for example _Fau vet, L. et al.| 
(|20TT|. 
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Figure 2. Polarisation fraction (top) and angle (bottom) in Galactic co-ordinates for our model of thermal dust emission at 150 GHz 
for the BSS (left column) and LSA (right column). The polarisation angle colour bar ranges from — 7r/2 to 7r/2. Both models include 
large (Is) and small (ss) scale magnetic field components. The ss turbulent component was added in the one dimensional, line-of-sight 
approximation and can be seen as an uncorrelated noise addition to the coherent Is component. There are significant differences in the 
morphology of the polarisation fraction between the BSS and LSA models due to the BSS model including the spiral arm structure. 



2.4 Stokes Parameters 

We combine the small-scale (ss) and large-scale (Is) magnetic 
field values according to 



Br 

Be 

B(h - 



- J^r,ss I -t>r,ls ; 

: Be^ss + Be^is , 

B(f)^SS + Bcf),ls , 



(4) 



where r, ^, and are now Solar-centric spherical polar co- 
ordinates. The polarisation at each point along the line-of- 
sight r is determined by the perpendicular field components, 
Be and 5^. 

The Stokes parameters for this model are projected 
out from our three-dimensional model using the appropri- 
ate line-of-sight integrals, 



Urr. 



nr 

Jo 

nr 
^0 



nd(r)dr , 

nd{r)po[B^{rf - Be{rf]dr , 

nd{r)pQ[2B^{r)Be{r)]dr , (5) 



and the normalisation po is set to reproduce the average 
polarisation fraction reported by WMAP outside their P06 
mask, 3.6% (Kogut et al. |2007 ). Here, e is the emissivity of 
the dust as a function of frequency, v. Note that we conform 



to the default convention applied in the healpi:?|j package 
( Gorski et al.||2Q05) regarding the sign of U . 

We have chosen the 3.6% average polarisation fraction 
as a reference value but the templates can be scaled to fit 
any other preferred value based on more detailed knowledge 
of the polarisation fraction in smaller patches of the sky. 
It is also useful to note that since we rescale the Q and U 
components the overall normalisation of the magnetic field 
model becomes irrelevant. However, the relative contribu- 
tions from the Is and ss components in the field remains as 
a model parameter. 

For the line-of-sight integrals we integrate from zero out 
to a maximum line-of-sight distance rmax of 30,000 pc. The 
integrals are discretised in steps of 0.1 pc. The direction of 
the lines-of-sight are chosen to coincide with the centre of 
all HEALPix pixels at a given A^g^de? where A^g^de is less than 
or equal to A^side of the total intensity template FDS map. 

From this model we require maps of the polarisation 
direction, 7, and degree, P, which are given by 



7(^,0) 



{Ql 



+ U; 



model/ 



arctan 



Jmodel 

^model 



(6) 



See http://healpix.jpl.nasa.gov 
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Figure 3. Stokes parameter maps (from top to bottom /, Q and U) in Galactic co-ordinates for our model of thermal dust emission at 
150 GHz for the BSS (left column) and LSA (right column). 



Figure [2] shows maps of P and 7 obtained from a line- 
of-sight integration at resolution A^g^de — 128 for the BSS 
and LSA magnetic field models including a one dimensional, 
small-scale turbulent component. The turbulent component 
is seen here as an uncorrelated noise contribution to the 
large-scale correlations induced by the large-scale magnetic 
field model. These maps can be compared to the "geometric 
suppression" factor shown in the right panel of Figure 8 of 
[Page et al.| (|2007| . There are significant differences between 
the BSS and LSA field models in the morphology of the 
polarisation fraction on the sky. The difference is greatest 
towards the Galactic centre and bulge and the Galactic anti- 
centre which coincides with a spiral arm. The LSA model 
does not include any azimuthal dependence and as such does 
not model any modulation of the magnetic field strength 
between spiral arms. In addition, the pitch angle of the LSA 



model, as fit to the WMAP data, is very low and this leads 
to a very mild dependence of the field alignment in the radial 
direction. These differences lead to a significantly simpler 
polarisation structure in the LSA model than in the BSS 
case which models the spiral arm structure explicitly. 

The final dust model at frequency v can be written as 



/du.t(e,</')=/FDs(e,^), 



(7) 



where /fds is the total intensity FDS prediction. 

Our final product is a template foreground map with 
small-scale structure modeled by the FDS predictions in the 
total intensity but with polarisation fraction and angle de- 
termined internally by our magnetic field model and line-of- 
sight integrals. An alternative approach taken by [Page et al.| 
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3D Polarization Amplitude, nside=1024 



ID Polarization Amplitude, nside=128 




0.0 ^^^ ^^ 4.0 /^K 
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Figure 4. Gnomic projection of the polarisation amplitude and angle in a 75 x 75 degree patch in the southern Galactic hemisphere. 
The left column shows the amplitude and angle from a line-of-sight integration including a full three dimensional realisation of the 
small-scale turbulent magnetic field model. A^J^^ = 1024 was used to calculate the "3D" maps but only for lines-of-sight corresponding 
to pixels inside the patch. The right column shows the same area from the full-sky templates with A^^^^ = 128. The full-sky maps used 
a one dimensional realisation of the turbulent component along the line-of-sight to speed up the computation. The absence of correlated 
small-scale structure and lower angular resolution of polarisation information in the "ID" case is clearly seen when comparing maps. 



( 2007 ) is to replace 7 with a map 7dust = 7* + 7r/2 where 7* 



is a smoothed map of observed starlight polarisation direc- 
tions. This approach, however, is limited by the resolution 
of the starlight data with only 1578 observations scattered 
around the sky. It also requires a large smoothing kernel of 
approximately 10 degrees in size and limits the application 
of any template derived in this way to very large scales on 
the sky, corresponding to angular multipoles ^ < 15, and 
Galactic latitudes \h\ > 10°. 



3 SCALES 

It is important to consider the range in angular scales our 
model is valid for. All our maps are pixelised at Aside = 1024, 
this ensures that the small-scale structure in the FDS pre- 
diction is oversampled since the IRAS resolution translates 
into a limit in angular multipoles of roughly ^fds "^ 1700 
and the healpix pixel smoothing scale is ^pix ^ 4 Aside- The 
overall, effective resolution of our templates is therefore lim- 
ited by the angular resolution of our line-of-sight coverage 
which is set by the healpix resolution Ag^^^. 

For the full-sky maps presented here and made available 
publicly we have chosen Ag^^e = 128 which corresponds to 
a limit of roughly i ^ 500 in multipole space. We also show 
in our example maps a small patch prediction with Ag^^e — 
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1024 (see Section p| which again oversamples the resolution 
given by the FDS templates. 

It is also important to consider physically relevant scales 
that enable the interpretation of the structure in our tem- 
plates. The most important of these is the injection scale 
for the turbulent, small-scale component of the magnetic 
field. We have set this to 100 pc. To obtain a rough estimate 
of the angular scales at which this physical scale becomes 
important we can use a "dust weighted" distance measure 
(r) = f r nd{T) dr / f nd{r) dr ~ 7000 pc for a mid-Galactic 
latitude line-of-sight. This can be used to place the angular 
multipole scale of injection at iinj ^^ 220 or roughly 1 degree. 
Beyond these scales the stochastic, turbulent component be- 
gins to dominate the structure in the polarisation and the 
model is only a statistical description of the real sky on 
these scales. An exhaustive exploration of foreground effects 
on scales below a degree would therefore require a Monte 
Carlo approach. 



4 MAPS 

We show a selection of template, full-sky maps at 150 GHz 
in Figure [3] /, Q, and U Stokes parameters are shown for 
both BSS and LSA derived templates (other frequencies are 
available on the on-line repository). The maps are of ther- 
modynamic CMB temperature in /jK units and are shown 
in Galactic co-ordinateq3 

As detailed above, the Q and U components have been 
normalised such that the average polarisation fraction out- 
side the area defined by the WMAP P06 mask is 3.6%. The 
resolution of the HEALPIX maps is A^side = 1024 but the po- 
larisation information is based on a line-of-sight integral at 
an angular resolution of A^s^de = 128. 

The maps have been obtained by the line-of-sight inte- 
gration of a magnetic field model that includes a small-scale 
turbulent realisation only along the line-of-sight direction, 
ie. our "one dimensional" approximation. Whilst computa- 
tionally intensive, "3D", full-sky maps that include a full 
three dimensional realisation of the turbulent component 
can be obtained, if required, with computation times of the 
order of 10 days. However we show results for a smaller 
75 X 75 degree patch in the southern Galactic hemisphere 
in Figure [4] These maps were obtained using a full three di- 
mensional realisation at an angular line-of-sight resolution of 
N^ide — 1024 and are compared with the same patch in the 
full-sky "ID" maps. The difference between the two is most 
clearly seen in comparing the polarisation angle which is un- 
correlated with the FDS intensity template. The full three 
dimensional case contains correlated structure on smaller 
scales due to the coloured power spectrum of the realisation. 
In contrast the one dimensional case is uncorrelated on small 
scales whilst preserving the large-scale correlations induced 
by the fixed, large-scale magnetic field model. Tailored, high- 
resolution, "3D" realisations of small patches such as those 
shown in Figure |4] are most useful for sub-orbital experi- 
ments that can only observe a limited fraction of the sky. 

^ Care must be taken in rotating Stokes parameters into other 
co-ordinate systems such as ecliptic and we have provided ro- 
tated maps on the on-line repository since most applications will 
simulate observations in this frame. 



5 CONCLUSIONS 

We have described a model of the polarised foreground that 
we expect to observe due to emission from dust within our 
Galaxy. The model uses a three-dimensional model of the 
Galactic magnetic field and dust field and integrates along 
the line-of-sight to each healpix pixel to obtain a polari- 
sation amplitude and angle. This information is combined 
with total intensity, FDS derived template maps at different 
frequencies to obtain a complete, polarisation template of 
foreground emission by interstellar dust. 

We have concentrated on two popular models for the 
structure of the large-scale structure of the Galactic mag- 
netic field, namely, the BSS and LSA models. The param- 
eters for the BSS model have been calibrated directly from 
measurements of the strength of the Galactic magnetic field. 
In the LSA case we have employed the parameters obtained 



by Page et al. (2007) in fitting to the WMAP observations. 



We calculate the polarisation alignment internally to our 
model in both cases since there is not sufficient external 
information on polarisation angles at resolutions relevant 
in this work. Some differences exist between the BSS and 
LSA derived templates but these are mostly at low Galac- 
tic latitudes away from the Galactic centre and as such ex- 
periments targeting small areas at high Galactic latitudes 
will not be sensitive to the differences. The differences do 
indicate however that a more accurate model of the Galac- 
tic magnetic field is required to produce realistic polarisa- 
tion templates for low Galactic latitudes. In the future, the 
Planck mission will provide an important test of Galac- 
tic magnetic field models through detailed characteristion 
of galactic foregrounds. 

We have developed a one dimensional approximation of 
the stochastic, turbulent, small-scale component of the field 
for obtaining full-sky templates. A full three dimensional 
realisation of the turbulent component can be used to obtain 
higher resolution templates for smaller patches of the sky. 

In future work we will be extending the model to in- 
clude synchrotron emission to form a complete picture of 
foreground emission relevant for polarisation experiments. 
Other developments will be required to increase the fidelity 
of the templates on small scales. These include the addition 
of a stochastic, small-scale density field to model small-scale 
structure in the density. In full "3D" calculations this will 
require the generation of an additional three dimensional, 
turbulence realisation which is correlated to the small-scale 
magnetic field. In addition, it would be useful to develop a 
simple model for the correlation of both realisations with 
the small-scale structure in the FDS derived total intensity 
templates. 

There is significant freedom in the parameters defining 
the small-scale structure in the templates. Experiments tar- 
geting small angular scales over small patches of the sky will 
be most sensitive to variations in the parameters and also 
to the stochasticity of the structure in the templates. Fur- 
ther Monte Carlo explorations of the variation in the maps 
is therefore warranted to quantify the impact of foregrounds 
on future sub-orbital experiments. As part of future work 
we will generate large ensembles of random realisations of 
the templates on small patches of the sky for the purpose of 
Monte Carlo studies. 
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The maps obtained from this model are available for 
use and can be downloaded from an on-line repositor}jj 
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